Superdense cosmological dark matter clumps 
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The formation and evolution of superdense clumps (or subhalos) is studied. Such clumps of 
dark matter (DM) can be produced by many mechanisms, most notably by spiky features in the 
spectrum of inflationary perturbations and by cosmological phase transitions. Being produced very 
early during the radiation dominated epoch, superdense clumps evolve as isolated objects. They 
do not belong to hierarchical structures for a long time after production, and therefore they are 
not destroyed by tidal interactions during the formation of larger structures. For DM particles 
with masses close to the electroweak (EW) mass scale, superdense clumps evolve towards a power- 
law density profile p(r) oc r~ 1,8 with a central core. Superdense clumps cannot be composed of 
standard neutralinos, since their annihilations would overproduce the diffuse gamma radiation. If 
the clumps are constituted of superheavy DM particles and develop a sufficiently large central 
density, the evolution of their central part can lead to a 'gravithermal catastrophe.' In such a case, 
the initial density profile turns into an isothermal profile with p oc r~ 2 and a new, much smaller 
core in the center. Superdense clumps can be observed by gamma radiation from DM annihilations 
and by gravitational wave detectors, while the production of primordial black holes and cascade 
nucleosynthesis constrain this scenario. 

PACS numbers: 12.60.Jv, 95.35.+d, 95.85.Pw, 98.35. Gi 



I. INTRODUCTION 



Gravitationally bound structures in the universe have 
developed from primordial density fluctuations S(x, t) = 
Sp/p that in turn were produced at inflation from quan- 
tum fluctuations. In the standard approach to infla- 
tion, the spectrum of these primordial fluctuations has 
a nearly scale-invariant form, P(k) = b\ oc k Up with 
n p ~ 1. During the radiation-dominated (RD) epoch 
fluctuations grow slowly, Sk oc ln(t/ti), while they grow 
as 5k oc (t/t eq ) 2 / 3 after the transition to the matter- 
dominated (MD) epoch at t = t cq . Gravitationally bound 
objects are formed and detach from the cosmological ex- 
pansion, when fluctuations enter the non-linear regime 
6 > 1. The non-linear stage of fluctuation growth has 
been studied both analytically [lJ-Ql and in numerical 
simulations 043 for the formation of galaxies and struc- 
tures on larger scales. The density profile in the inner 
part of dark matter (DM) halos is given by p(r) oc r^ 13 , 
with (3 w 1.7 — 1.9 in analytic calculations [3], /? = 1 
in the simulations of Navarro, Frank and White [fj| and 
/3 = 1.5 in the simulations of Moore et al. [|[ and Jing 
and Suto [<|. 

The smallest DM objects in the universe, which we 
shall call clumps or subhalos, are produced first. The 
evolution of DM clumps has been studied in Ref. [10] in 
the hierarchical model in which due to the merging of ob- 
jects a small clump is hosted by a bigger one, the latter 
is submerged into an even bigger one, etc. The impor- 
tant observation of [10( was the role of tidal interactions, 
which fully disrupt most clumps. The survived clumps 
can be further destroyed in the Galaxy by tidal interac- 



tions in the Galactic plane, near the Galactic center, and 
in collisions with stars in the halo (see [Tl| for a review). 
The characteristic feature of these processes of disrup- 
tion is that the core of a clump survives and thus the 
gamma signal from DM annihilations in clumps changes 
only mildly [ll|. A statistical approach to the search 
for galactic small-scale substructures has been recently 
proposed in [l2l lllij. 

The mass spectrum of DM clumps has a low-mass cut- 
off M m i n due to the leakage of particles from a clump. 
This mass is strongly model dependent: It depends on 
the leakage mechanism (free streaming, collisional damp- 
ing, etc.), on the properties of the DM particles and 
the resulting decoupling temperature and others. There- 
fore, the predicted M m i n varies for neutralinos in the 
minimal supersymmetric standard model from 10 -7 to 

io- 5 M Q Mm- 

We have described above the standard cosmological 
scenario for the clumps. In non-standard scenarios the 
properties of DM clumps can be very different. In 
Ref. [16j . isothermal perturbations in the DM density 
were considered within the framework of a spherical col- 
lapse model. Perturbations collapse in the RD epoch and 
produce superdense DM objects. Another possibility for 
the production of superdense clumps is given by a spiky 
spectrum of perturbations [i~^ - fl9l ]. The general idea com- 
mon to these scenarios is that there exists a spike on top 
of a scale-invariant power-law spectrum of perturbations 
which results in the production of dense clumps in a very 
early cosmological epoch. In this work, we consider in 
contrast to [16| the formation of clumps at the RD epoch 
from adiabatic spiky perturbations. The difference to 
isothermal perturbations is mainly in their evolution dur- 
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ing the linear stage: While isothermal perturbations are 
frozen in, adiabatic fluctuations grow logarithmically. 

We include in this and in the accompanying paper [20I ] 
a discussion of the detection prospects of stable super- 
heavy DM particles. Since the annihilation signal from 
the mean distribution of these particles in the halo is 
far below observational limits, we examine whether there 
are new effects which improve the detection chances. One 
such effect follows from the early kinetic decoupling of su- 
perheavy DM particles from the thermal plasma. In this 
case the cutoff mass can be significantly smaller, as e.g. 
in the case of ultra-cold WIMPs [l8| , and thus clumps of 
practically all masses are formed. This opens the door 
for the formation of light superdense DM clumps at the 
RD stage. The only necessary condition is the existence 
of spiky small-scale perturbations. 

This article is organized as follows. We determine the 
initial properties of the DM clumps, first assuming a 
standard power-law for the initial cosmological pertur- 
bations in Sec. |TT] and then a spiky perturbation spec- 
trum in Sec. IIIII In Sec. IIV1 we derive constraints on the 
superdense clump scenario considering primordial black 
hole production. Then we study the evolution of the 
density profile of superdense clump in Sec. lVIi comment- 
ing on the case of neutralinos with masses close to the 
electroweak mass scale in Sec. |Vj We present finally our 
conclusions in Sec. IVIII 



II. CLUMPS IN THE STANDARD 
COSMOLOGICAL SCENARIO 

We briefly remind in this section the formation of 
clumps and their properties assuming a standard power- 
law spectrum of the initial cosmological perturbations. 
In contrast to the usual approach, we allow here very 
small masses of the clumps being inspired by the small- 
ness of M m i n in the case of superheavy DM (SHDM), 
where M m ; n can be of order of SHDM particle mass m. 

Small clumps form at the MD epoch for z ^> 1, i. e. at 
a time when the effect of the cosmological constant can 
be still neglected. In the spherical model of the Press- 
Schechter theory [2l|, [22j], the formation of an object oc- 
curs at the time tf when the density contrast 5(M,tf) 
reaches S c = 3(12tt) 2 / 3 /20 ~ 1.686. The mean density 
Pint and the radius R of the collapsing clumps are 



Pint = Kp(Zf) = Kp cq ( -!- L j = K p c 
\ t + Z cq J 
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R 



3M 

47T / 5inf 



1/3 



(1) 

(2) 




log M, g 



FIG. 1: The mean density p of DM clumps as function of 
the clumps mass M for different spectral indices n p of the 
primordial density perturbations. 
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FIG. 2: The fraction of DM in the form of survived clumps 
per unit logarithmic mass interval SM ~ M as function of 
of clump mass M for n p = 0.949 (bottom), 0.963, 0.977 and 
1 (top): The initial fractions are shown by dashed lines, the 
present fractions by solid lines. 



According to Ref. [11], surviving clumps are charac- 
terized by v ~ 1 — 3 and we set v — 2 in all following 
calculations. Having fixed v, the dependencies R(M) and 
p(M) are unambiguous and the mean density p of small- 
scale DM clumps as function of the clumps mass M is 
shown in Fig. [T] 

The mass function of clumps, i.e. the fraction of DM 
in the form of clumps with mass M, is given by fllj ] 



dM 
tint M 



0.02(n 



dM 
6) M ' 



(3) 



where k = 18n 2 ~ 178 [22|], a eq (M) is the variance and 
v = 5 eq / a eq (M) is the peak height of the density fluctu- 
ations at the time t eq of matter and radiation equality, k np with 



where the effective exponent n in Eq. ([3]) is found as 
n = —3(1 — 2<91og(T cq (Af)/<91ogM) and depends very 
weakly on M . The simplest inflation models give P(k) ex 



while p C q is the density at t e 



w 1. The 7-year WMAP data, n p = 
0.963 ± 0.014, favor clearly n s < 1 23]. Clumps can 
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form nevertheless, because of the presence of additional 
logarithms in the transfer function. The small-scale spec- 
trum at the epoch of matter-radiation equality can be 
expressed as 



cpo 
10] 



cr eq (M) ~ 8.2 x io 3 - 7 ^- 1 )- 3 



M 



Mr. 



Q. 
CM 



(4) 



The mass function ([3]) with the spectrum Q is shown 
in Fig.[5]by dashed lines. Its 1/M shape is in good agree- 
ment with the corresponding numerical simulations of 
Ref. [24j |. only its normalization is a few times smaller 
than the one found there. For an extrapolation by many 
orders of magnitudes this must be considered as remark- 
able agreement. 

Note also that using the power-law spectrum that is 
normalized to the temperature fluctuations of the CMB, 
i.e. at cosmological scales, for sub-galactic scales or even 
DM clumps with mass M ~ lg implies an extrapolation 
by ~ 48 orders of magnitudes. This extrapolation can be 
justified only within the simplest models for inflation. 

Integrating the mass function ([3]) from M mm to M ~ 
10 2 M Q , we obtain the initial (i.e. before possible destruc- 
tion in the Galaxy) fraction of DM in the form of clumps. 
In contrast to the standard case of EW scale neutrali- 
nos, where M min ~(KT 6 -1CT 8 )M Q (lSf. in superdense 
clumps the DM particle can have much larger mass and 
thus M m i n can be much smaller. As a result the frac- 
tion of surviving clumps increases. In particular for su- 
perheavy neutralinos, A/ m i n can be comparable to the 
particle mass m and the fraction of surviving clumps is 
calculated as as 0.15, 0.18, 0.21 and 0.26 for n p = 0.949, 
0.963, 0.977 and 1, respectively. 

Clumps inside galaxies lose mass and can be destroyed 
in tidal interactions with stars. The collective gravita- 
tional field of the Galactic disk is the most important 
factor for the clump destruction. A method to study the 
destruction process of clumps was presented in fill ] (for a 
more detailed approach with gradual mass loss see }lll|). 
where only clumps with M > 10 -6 M Q were considered. 
Here we calculate the survival probability for the wider 
mass interval m < M < 10 2 M Q , using the same formal- 
ism as in [llj]. The result for the survival probability P(p) 
at the position of the Sun, r — 8.5 kpc from the Galac- 
tic center is presented in Fig. [3j Note that the survival 
probability P(p) means the fraction of surviving clumps 
near the Sun but most of these clumps have elongated 
orbits and spend the largest part of their orbital period 
far from the Sun at the outer parts of the Galactic halo. 

The resulting mass function P£i n t that accounts for the 
effect of tidal destruction by stars is shown in Fig. [2] by 
solid lines. Integrating J P£i nt dM/M again from M ~ m 
to M ~ 10 2 M Q we obtain the actual fractions of DM in 
the form of clumps as 0.006, 0.015, 0.033 and 0.085 for 
n p = 0.949, 0.963, 0.977 and 1, respectively. 
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FIG. 3: The survival probability P(p) as function of the mean 
internal clump density p at the distance 8.5 kpc from Galactic 
center. 



Clumps formed from the standard power-law spectrum 
considered above have a rather small density. For many 
DM particle candidates, including SHDM particles, such 
clumps are unobservable via their annihilation signal and 
these clumps can be detected only gravitationally. It 
has been already suggested that interferometric detec- 
tors for gravitational waves like LISA have the capa- 
bility to detect the tiny variation of the gravitational 
field, when a compact object crosses the detector. Small 
SHDM clumps should be included in the list of objects 
to be searched for by LISA, such as primordial black 
holes [25|], asteroids [26| or compact DM objects of un- 
known nature (27j . The observable signal is caused by 
the gravitational tidal force which changes the interfer- 
ometer arm length and produces correspondingly a phase 
shift. LISA will have the capability to search for com- 
pact objects in the mass interval 10 16 g< M < 10 20 g 
according to Ref. [2j| and 10 14 g< M < 10 20 g accord- 
ing to Ref. [27| ■ The signal will be in the form of sin- 
gle pulses with its characteristic frequency at the lower 
end of the expected LISA sensitivity curve and a rate 
~ a few per decade, if the objects constitute the ma- 
jor part of DM. The clumps under consideration in this 
Section present only 1-10% of all DM, and correspond- 
ingly, the detection rate will be 1-2 order of magnitudes 
lower. In addition, the radii of the clumps generally ex- 
ceed LISA's arm length L ~ 5 • 10 11 cm (see the Fig. 2]) 
and the tidal forces will be smaller due to the extension 
of these objects. Therefore, the detection of the SHDM 
clumps by LISA seems unlikely. The next generation of 
gravitational wave interferometers offers more promising 
perspectives for detection (for details see Ref. [25[). 
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FIG. 4: The mean virial radius R of DM clumps ([2]) as func- 
tion of the clump mass M for several values of the spectral 
index of primordial density perturbations (from top to bot- 
tom): n p = 0.949, 0.963, 0.977 and 1. The horizontal line 
shows LISA's arm's length L ~ 5 • 10 11 cm. 



proposed also in [34J . The authors of Ref. [35J found ev- 
idence for excess power at small scales ~ 10/i -1 kpc in 
comparison with a flat primordial power spectrum. This 
result was obtained from the study of Lyman-a absorbers 
and can be explained within complex inflation models 
with the generation of extra power at small scales. Such 
models can lead to the effective production of very dense 
clumps. 

We will refer to all these models collectively as spiky 
models or spiky mass-spectrum models. 

Dark matter clumps are formed in a wide range of 
masses, if the power spectrum of primordial cosmolog- 
ical density perturbations has a power-law form. If on 
the contrary the spectrum has a peak on some scale, then 
clumps are formed mostly in a narrow range of masses, 
near the mass that corresponds to this peak. 



B. Formation of superdense DM clumps at the RD 

epoch 



III. NON-STANDARD PERTURBATIONS AND 
SUPERDENSE CLUMPS 

A. Spiky density perturbation spectrum 

The variance of the normalized power-law spectrum at 
the horizon scale during the RD stage was expressed for 
the standard inflationary scenario in Ref. [28| as 



a H {M) ~ 9.5 x 10" 



M 



10 56 g 



(5) 



We see that in view of current observations (n p < 1) the 
variance an {M ) is too small for the formation of clumps 
at the RD stage. Such clumps can be produced effectively 
only from non-standard spectra containing e.g. spikes. 

A sharp peak emerges in the fluctuation spectrum 
e.g. , if an inflationary potential V(tf>) has a flat segment 
[l7l [29l | . The mean density perturbation on the hori- 
zon scale is <5h ~ Mp 3 V 3 ^ 2 /V . Hence, if the derivative 
V = dV \<j>) I d(j) — > for some value of the scalar field 
<f>, then a peak emerges in the perturbation spectrum 
on the corresponding scale. A similar effect can arise 
in inflationary models with several scalar fields [3(| [HJ ■ 
In both types of models, the spectrum outside the peak 
can have an ordinary shape. In particular, it can be a 
Harrison-Zel'dovich spectrum, and can give rise to galax- 
ies, clusters and superclusters according to the standard 
scenario. 

Another possibility to generate a spiky density pertur- 
bation spectrum are cosmological phase transitions, for 
example the QCD phase transition 32]. If somewhere 
a high peak arises in the perturbation spectrum, then 
the corresponding clumps would be the densest DM ob- 
jects in the universe. Theoretical models for nonstandard 
spectra were discussed also in (33[. A peak in P{k) was 



A useful approximation for the nonlinear evolution of 
perturbations in the radiation dominated epoch is the 
spherical collapse model 0, [36j]. In this model, the evo- 
lution of perturbations after the horizon crossing is de- 
scribed by 



, 3 

i + - 2 y 



db 

dy 



1 + $ 



= 0, (6) 



where y — a(ry)/a cq , r\ = dt/da is the conformal time, 
a cq is the scale factor at rj cq , and $ = ^pdm/pdm is the 
relative overdensity of DM. The radius of the perturbed 
region is parametrized as 



r = a(r))b(r))£ . 



(7) 



Here, £ is the comoving coordinate of the spherical layer 
considered and the value b(rj) takes into account the slow- 
down of the cosmological expansion in the perturbed den- 
sity region. Equation (|6]) is applicable for the evolution 
of both entropy and adiabatic perturbations, but has to 
be used with different initial conditions. 

The formation of clumps from entropy perturbations 
was considered in fl6| . In this particular case, the initial 
data have the form $ = 5pdm / Pt)m and db/dt = 0. The 
object formed has the density [n| 



p~ 140$ 3 ($ + l)p c 



(8) 



For instance, $ ~ 1-^10 4 in the case of axions as DM, and 
axionic miniclusters have masses in the range ~ (10~ 13 -f- 
0.1)M Q . The observational signatures of the presence 
of these axionic miniclusters in the Galactic halo were 
considered in [l6|, [3?| • 

The corresponding method for the nonlinear evolution 
of adiabatic perturbations during the radiation domi- 
nated epoch is described in 38]. For adiabatic pertur- 
bations $ = 0, the initial velocity db/dt is nonzero and 
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is defined using linear perturbation theory. The trans- 
formation from the Euler description for the growth of 
density perturbations 5 to the Lagrange description (0 
is provided by the relation b = (1 + 5)" 1 / 3 [36]. The 
evolution of perturbations with 5 <C 1 on scales less than 
the horizon is defined by the known analytic solution [32[ 
(see also (38J) 



6 = 



3A m 



IE 



(9) 



In this solution the numerical constant equals 7^ — 1/2 
0.077, A in — 5ft/(f), (f) — 0.817, (5 H is the radiation density 
perturbation on the horizon scale and the variable x is 
related to the comoving wave- vector k of the perturbation 
by x — krj. The connection between x and y is defined 
by the relation [38| 
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M 1 / 3 G 1 / 2 pli 6 



(10) 



It is suitable to connect the analytic solution of the lin- 
ear theory ([9]) with the numerical solution of the nonlin- 
ear Eq. <j6j) at the time corresponding to the "transition" 
value of perturbations with 8 — 0.2 (see [38]). At this 
moment we define the initial velocity of the forming DM 
clump as 



db 
dy 



Sub 4 
" 2 y( f> 



(11) 



The cosmological expansion of the forming DM clump 
stops when dr/dt — or according to Eq. ([7]) when 
db/dy = —b/y. The corresponding density and radius 
of the clump are 



Pmax — Peq2/ ma x^maxi Rmax 



/ 3M 



1/3 



(12) 



where 6 ma x and y ma ,x are respectively the values of b and 
y at the same moment. After decoupling from the cosmo- 
logical expansion, the object virializes and contracts by 
a factor two. In Ref. [38] this model was used to describe 
a noncompact DM object with single mass ~ 0.1M Q , 
presumably observable through microlensing. Now we 
consider the whole possible range of masses and densi- 
ties of DM clumps. Calculating numerically the solution 
of Eq. ([6]) within the above formalism, we find the den- 
sity of the clump p — p(M,5u) as function of its mass 
M and the radiation perturbation value on the horizon 
scale 5n as shown in Fig. [SJ 

Some characteristic values of the clump density p 
are displayed in Fig. [6] for several values of the clump 
mass M. One observes the convergence of curves to 
p ~ p cq ~ 10~ 19 g cm~ 3 at small <5h, i-e. for clumps 
formed near matter-radiation equality. This corresponds 
to the known analytical results that the evolution during 
the MD epoch does not depend on the mass but only on 
the initial (at t = t eq ) value of the fluctuation. 
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FIG. 5: The mean density p (in g cm" 3 ) of DM clumps as 
function of the perturbation dn in the radiation density on 
the horizon scale and the clump mass M (in M@). 
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FIG. 6: The mean density p (in g cm -3 ) of DM clumps as 
function of the perturbation Sn in the radiation density on 
the horizon scale; solid lines from top to bottom are for clump 
masses M = 10" 10 , 1, 10 10 , 10 20 , 10 30 g. The dashed line is 
the bound on the clump density from primordial black holes 
overproduction with threshold 5 C = 0.7. The time of two- 
body gravitational relaxation inside the clump core is less 
than the age of the Universe for clumps above the dotted 
lines, if the DM particle mass is m > 10 11 GeV. 



Note that (in contrast to the case with standard power- 
law spectrum of cosmological perturbations) superdense 
clumps from a spike in the spectrum are not destroyed 
by tidal forces and their mass function peaks near a def- 
inite mass. Therefore the fraction of DM in the form 
of such clumps is £ ~ 1/2. Half of the volume is 
in the form of overdensities (clumps), and the remain- 
ing space is filled by voids. Because of the compact- 
ness of superdense clumps, these clumps can satisfy the 
condition R < L ~ 5-10 11 cm for the mass interval 
10 g< M < 10 20 g and are thus observable by the 
LISA detector. 



IV. CLUMPS AND PRIMORDIAL BLACK 
HOLES 



The fraction of the mass in radiation that is trans- 
formed into PBHs at the time in is then [4(| 



The formation of DM clumps leads to several restric- 
tions on the fluctuation spectrum. For instance, high- 
energy particles from DM annihilations in clumps during 
the epoch of nucleosynthesis and after it (the cascade 
nucleosynthesis) might distort the prediction of standard 
nucleosynthesis . 

Another important restriction on the spectrum of 
the adiabatic perturbations comes from upper limits on 
the mass and density of primordial black holes (PBHs) 
[39l lioj , because the value of DM density perturbations 
depends on the radiation density perturbations and the 
formation of DM clumps can be related to the formation 
of PBH from the same perturbation spectrum 38[. It 



should be noted that in the case of entropy perturbations 
PBHs do not form. 

Clumps and PBHs originate from fluctuations of the 
same type but are formed at different times. The large 
difference between the masses of DM clumps and of PBH 
arises because of the large difference in energy density 
enclosed in a fixed comoving volume as function of time: 
The energy density of radiation at the RD epoch far ex- 
ceeds the mass in CDM at the matter domination epoch 
in the same comoving volume. 

The formation of PBHs takes place on the tail of the 
distribution of Gaussian fluctuations, whereas the main 
part of clumps is produced from r.m.s. fluctuations. 
Therefore only a small part of the fluctuations which re- 
sult in the formation of clumps may produce PBHs at the 
RD epoch. In other words, because of the large thresh- 
old of PBH formation, the major part of fluctuations 
does not collapse into PBHs and evolves continuously 
up to the end of the RD epoch. During the RD epoch 
the mass of radiation in the comoving volume varies as 
M r (t) = Ma(t cq ) /a(t), where the scale- factor of the Uni- 
verse a oc t 1 / 2 and M is the comoving mass at the mo- 
ment of transition to MD. The mass M equals approxi- 
mately to the mass of non-relativistic matter inside the 
fluctuation, i.e. the mass of a clump which may be formed 
from this fluctuation. On the other hand at horizon cross- 
ing M r (t) ~ 47r(2ci) 3 ( o(i)/3, where pit) = 3/32nGt 2 . 
^From these relations for M r (t) we estimate the mass 
Mh and the time tn of PBH formation as function of 
the clump mass M as 
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H 



cM 2 ' 3 G- 1/3 t 1 J 3 t H ~ GM H /c 3 



(13) 



From the Friedmann equations, the formula for Mh was 
derived exactly in Ref. [38] and is given by 



W 3 



1/6 



M 2 ' 3 c 



while the dependence tu(M) is 



= 2x10 s 
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M r . 



(14) 
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where S c is the threshold value of the density perturba- 
tions Sh which result in PBHs formation. The current 
PBHs density parameter is J7bh — ^o,(t cq )/a{tn). 

For a large enough value of the r.m.s. perturbation 
Ah = (^h) 1 ^ 2 ' an extremely large number of PBHs can 
be formed 40] . This provides a limitation on Ah • 

The number density of PBHs depends strictl y o n the 
threshold value S c . In early works, e. g. [40l - l42| . the 
value of S c = 1/3 was obtained. In recent years the phe- 
nomenon of critical gravitational collapse was discovered 
in numerical simulations, for which S c ~ 0.7 (43l . l44j . 
Some limits on the number density of PBHs in different 
mass ranges were obtained in [40l . [l2[ . These restrictions 
on the value of Ah for PBHs are shown in Fig. 2 for the 
case S c = 0.7. The relation (TT4"]) was used in our calcu- 
lations. The local minimum on the curve corresponds to 
the restrictions on the Hawking evaporating PBHs with 
masses Mbh — 10 15 g. For PBHs with a larger mass 
the only restrictions comes from the condition that their 
cosmological density parameter Opbh < 1- 

We recall that PBHs are formed on the tail of the Gaus- 
sian perturbation distribution, <5h > S c S> Ah- On the 
contrary the overwhelming number of DM clumps are 
formed from the r.m.s perturbations. For this reason in 
Fig. 2 and in the calculations for DM clumps we put 
Sh - A H - 



V. SUPERDENSE CLUMPS FROM ORDINARY 
NEUTRALINOS 

We consider in this section the standard case of ther- 
mally produced neutralinos with mass close to the elec- 
troweak mass scale. We will show that the diffuse gamma 
flux produced by such neutralinos constituting super- 
dense clumps exceeds the observed flux, and thus su- 
perdense clumps should consist of DM particles non- 
thermally produced. 

In order to make our estimate most transparent, we 
consider first the integral photon flux produced by DM 
annihilations. This flux is easy to estimate using the an- 
nihilation cross section (cry) for the process x + X ~ * 
7r°+all and the mean density p- ult of neutralinos in a 
clump. 

We calculate first the rate N 7 of gamma-rays with en- 
ergies higher than 70 MeV produced by a single clump, 
assuming a r~ 1,8 density profile with core at r < R c for 
a clump with total mass M and radius R, 
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where r/ w o is neutral pion multiplicity, R c = x c R is the 
core radius, and m x is the neutralino mass. 

The total diffuse flux produced in a galactic DM halo 
can be calculated as 



JZ 



■^-n c \R h N 1 



47T 



(18) 



with Rh as the radius of the DM halo and n c \ as the 
clump mean space density, which is given by the fraction 
£ of DM in the form of clumps and the mass Mh of the 
galactic DM halo as 



n c i 



3 jM h 1 
4tt M R\ 



(19) 



Using Eqs. (|17M19p we can express the diffuse flux J* ot in 
terms of the mean density of neutralinos in clumps p; nt , 
the main characteristic of superdense clumps: 



/nfw 



0ATi^R h p int p^f o 



(20) 



where pE*f Q is the mean density of DM halo. To take into 
account the NFW density profile one must multiply the 
homogeneous halo result by the additional factor /nfw = 
293. The obtained flux is given as convenient expression 
where most of parameters are observationally known and 
the main characteristic of superheavy clumps is p m t . For 
the parameters in Eq. (|20p . we use r\^o =10 appropriate 
for gauge boson decays, £ = 1/2 (for spiky scenario), 
R h « 200 kpc and pgfc = 1.1 x 10~ 3 GeV/cm 3 obtained 
as ZM h j ^R\. We assume i c »01 

Typically superdense clumps have very large densities 
(see Fig. [BJ and for ordinary neutralinos the resulting 
gamma-ray flux exceeds the observations. First we ana- 
lyze the problem, whether ordinary neutralinos are com- 
patible with a spiky scenario of clump production. With 
this aim we choose in Eq. (1201) parameters which min- 
imize the flux. For the mean density of neutralinos in 
a clump, pi n t, we take the minimum value, assuming 
neutralino produced at the beginning of the MD epoch. 
In this case p- mt = 178p eq = 8.3 x 10 6 GeV/cm 3 (see 
Eq. ((T|) in the limit Zf — » z eq ). We parametrize the an- 
nihilation cross-section (av) by the characteristic value 
1 x 10~ 26 cm 3 /s as 



M26 = M/(io- 



~ 26 cmV 



(21) 



With this parameters the minimum gamma-ray flux is 



Jl 



4.3(<to) 2 6TO 10 2 



(22) 



where toioo is the neutralino mass m x in units of 100 
GeV. 

The integral flux (|22|) is about 5 orders of magnitudes 
larger than the observed flux. Does it help to increase 
the neutralino mass or to consider a smaller annihilation 
cross section? To answer the first part of this question, 



we consider now the differential isotropic diffuse photon 
flux observed by Fermi-LAT at |6| > 60 degrees [50j . 



E \ 



-2.45 



-2 „-! gr -l 



J ohs (E) = 6 x 10- 7 I — j GeV" 1 cm- 2 s 

The differential photon flux produced by annihilations 
can be obtained from Eq. (f20| replacing 2r] 7r a by 
dN/(2m x dx), where dN/dx is the number of photons 
with energy E — xm x /2 produced per annihilation. 
Since moreover dN/dx increases at small x for increasing 
m x , the ratio J^(E)/J Q \, S (E) is practically constant. The 
minimal allowed annihilation cross-section of neutralinos 
obtained in [49| for the case of strongly suppressed s-wave 



annihilations is (av) — 1.7 x 10 



-30. 



"100 



cm 3 /s. With 



these parameters the minimum gamma-ray flux is still 
above the measured Fermi-LAT flux. Finally, we remark 
that even postulating at tree-level only couplings e.g. to 
electrons would lead to an overproduction of photons via 
Bremsstrahlung, c.f. e.g. Ref. [5l|]. In conclusion, ordi- 
nary neutralino as other thermally produced DM parti- 
cles are excluded as constituents of superdense clumps. 
A DM particle suitable to compose superdense clumps 
must have a smaller annihilation cross section than al- 
lowed for a thermal relic. 



VI. DENSITY PROFILE EVOLUTION 

In the case of a spiky spectrum, clumps are formed 
not in the process of hierarchical clustering but due to the 
evolution of isolated density fluctuations. Such a scenario 
is similar to the analytic approach in Ref. Q. The ordi- 
nary gravitational contraction combined with the multi- 
stream instability produces the universal power-law den- 
sity profile with exponent j3 = 1.7 — 1.9 [3j. This power- 
law shape for p(r) has been recently confirmed in the 
numerical simulations [45j for neutralino clump forma- 
tion during the MD epoch. We assume here that the 
clumps produced at the RD stage in the process of ordi- 
nary gravitational contraction have a profile p(r) oc r -1,8 
for R c < r < R and p(r) = p c = const for r < R c , where 
R c is the unknown core radius of the clump. This core 
may be produced due to tidal forces [l(| in the clumps 
formed at the MD stage. In this case a large core is pro- 
duced with x c = R c /R ~ 0.01 — 0.1. More precisely, the 
given value corresponds not to the radius of the constant- 
density core but to the break in the slope of the density 
profile. Moreover, the above-mentioned calculations are 
valid for the MD dominated epoch, where the process of 
core formation can be much different from that at the 
RD epoch. 

Another estimate for the core size has been obtained 
in Ref. [|[, where x c is defined by the damping mode of 
the perturbations. The authors obtained x c ~ <5 3 q , where 
S cq is the value of density fluctuation at the beginning of 
the MD stage. However, this estimate is also valid for 
the MD epoch. 
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At the current level of knowledge, the relative radius 
x c = R c /R of the core produced by ordinary gravita- 
tional contraction must be considered as a free parame- 
ter. In the most conservative case we use x c ~ 0.1. The 
central density p c depends on the mean clump density 
p = 8/9max (see (d2J) as p c = p/{Zx 2 c ). 

We shall briefly discuss the evolution of superdense 
clumps formed from superheavy particles. Quantita- 
tively, it will be considered in the accompanying paper 
II. 

The first stage of evolution, the ordinary gravitational 
contraction, proceeds like in the case of ordinary neu- 
tralinos and results in the production of a p(r) ~ p -1 - 8 
profile with a relatively large core, x c ~ 0.01 — 0.1. Other 
processes can become important at the second stage: (i) 
two-body gravitational scattering and (ii) some limit- 
ing effect like Fermi degeneracy or the intensive anni- 
hilation of particles. In the cores of superdense clumps 
with large densities n of particles the binary gravita- 
tional scattering of constituent DM particles with large 
masses m may become the dominant process, which 
causes the "gravithermal instability" or "gravithermal 
catastrophe" , well known in theory of globular star clus- 
ters. Note that this effect takes place only for super- 
heavy DM and only in the most dense parts of super- 
dense clumps. In Fig. [SJ this region is located above the 
dotted line. 

How can it be that gravitational two-body scattering 
becomes the dominant process? It occurs because gravi- 
tational scattering is proportional to m 2 , while EW scat- 
tering of these particles is proportional to 1/m 2 . The 
other two factors are the large density n of particles in 
the core and the long-range character of gravitational 
interactions. All this provides the fast gravitational re- 
laxation of the system. As a result of the gravithermal 
instability a clump develops an isothermal density profile 
p(r) cx r~ 2 with a tiny core. This core can be produced 
by the pressure of a degenerate gas in the case of super- 
heavy fermionic particles or by the inverse flow due to 
the annihilation of particles in the clump center (4|| Hi| . 
In these cases the radius of the new core is determined 
by the elementary particle properties of dark matter. 

VII. CONCLUSIONS 

Superdense clumps can be produced from isothermal 
perturbations [l6j or from spikes in the spectrum of adia- 
batic perturbations [13, IH[ . These objects are produced 
in the very early universe during the RD epoch. In princi- 



ple, the perturbation spectrum may include both a scale- 
invariant power-law component and spikes. Being pro- 
duced very early during the radiation dominated epoch, 
superdense clumps evolve as isolated objects. They do 
not belong to hierarchical structures for a long time after 
production, and therefore they are not destroyed by tidal 
interactions during the formation of large-scale struc- 
tures. 

In the case of EW scale mass particles, e.g. ordinary 
neutralinos, the density profile has a r -1 ' 8 shape with a 
relatively large core characterized by R c /R ~ 0.01 — 0.1, 
produced by tidal forces. Ordinary neutralinos are ex- 
cluded as the constituents of superdense clumps, because 
they overproduce the diffuse gamma-ray spectrum above 
100 MeV. The constituent DM particles in superdense 
clumps must be either very weakly annihilating or be su- 
perheavy, or both. The limit on the superdense clumps is 
imposed by primordial black holes which originated from 
the same perturbation spectrum. The allowed intrinsic 
densities of superdense clumps are shown in Fig. O The 
formation of superdense clumps at the RD epoch was 
studied previously using somewhat different assumptions 
in Refs. [II 11, [Hi. 

The density profile in superdense clumps depends on 
the properties of the DM particles. For very heavy 
constituent particles and large intrinsic densities of the 
clumps a "gravithermal catastrophe" (instability) may 
develop in superdense clumps. As a result the initial den- 
sity profile turns into an isothermal one, p- ln t (r) cx 1 jr 2 , 
and the large initial core collapses into a tiny, very dense 
new core. The steep density profile and the smallness of 
the core lead to a strong DM annihilation signal. The ra- 
diation produced by DM annihilations restricts this sce- 
nario, e.g. due to the cascade nucleosynthesis following 
standard nucleosynthesis. On the positive side, super- 
dense clumps can lead to detectable gamma radiation 
even in the case of superheavy DM particles [2(| ■ Super- 
dense clumps can be in principle observed also by gravi- 
tational wave detectors. 
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